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A new Fermionic Molecular Dynamics (FMD) model based on a Skyrme functional is 
proposed in this paper. After introducing the basic formalism, some first applications to 
nuclear structure and nuclear thermodynamics are presented. 



1. Introduction : the Fermionic Molecular Dynamics model 

The main goal of the Fermionic Molecular Dynamics model proposed by S. Drozdz et 
aL^and B. Caurier et a/P} and improved thereafter by H. Feldmeier^, is to modelise 
the dynamics of a fermions system where antisymetrisation is exactly taken into 
account. The FMD model consists of a single Slater determinant parametrised by 
a set of dynamical variational parameters Q(t) — [q^(t)\[i = 1,2,...] : 

io(*))=4nift(*)> a) 

k 

where the single particle states are chosen to be gaussians, in order to localise the 
particles in phase space and be able to deal with particle density fluctuations and 
clustering: 

:' r '//.<") = ^i> i - v 2ak[t) J I ui/ >.",,(/)) /»/(/.)} (2) 

The real parameters \k an d 4>k are the dynamical phases for the evolution of the 
spin degree of freedom, \m(k)) is the isospin degree of freedom, and the complex 
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parameters ak and bk are linked to the classical coordinates of nucleons. Indeed 
the expectation value of the position and the momentum for a single particle state 
reads : 

lm(a k (t)) 
He(o fc (t)) 



< r > = 



r kW = Ite (i W ) + ^W|i m (i(*)) (3) 
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The dynamical evolution of each variational parameters is obtained imposing a 
minimisation of the action : 

S f 2 dt{Q(t)\ih^ - H\Q(t)) = (5) 

where H is the Hamiltonian of the system. This equation gives the FMD equation 
of motion for each variational parameter : 

where TL = (Q(t)\H\Q(i)). An important feature of the FMD model is that the 
set of dynamical variational parameters includes the widths of the single par- 
ticle states. This means that the quantum-mechanical spreading of a free wave 
packet can be reproduced by the model. Another important point is that the 
single Slater approximation implies that the expectation value of the energy is 
given by a functional of the one body-density only, namely the Hartree-Fock energy 
(Q(t)\H\Q(t)) = E HF [f(p)}. The Hilbert space is restricted respect to TDHF be- 
cause of the gaussian wave packet constraint, meaning that FMD can be viewed as 
an approximation of TDHF. On the other side, it is interesting to remark that the 
FMD model gives the exact Hamilton equations for classical particles, meaning that 
all the correlations and fluctuations are taken into account at the classical level. 
Since the expectation value of the energy for the FMD model can be written as 
a density functional, the lastest Skyrme functionals as the SLy parametrisationfl 
can be used, with appropriate restoration of Galilean invariance 



2. FMD for nuclear structure 

Our FMD model based on a Skyrme functional can be tested comparing observables 
with HF calculations and experimental data. In the following we show calculations 
of ground state properties and Giant Monopole Resonance (GMR) frequencies. 
Concerning ground states, the FMD variational principle coincides with the mini- 
mization of the energy j^- = 0. Then FMD becomes a static model where = 0. 
To obtain the ground state properites, we have used an iterative algorithm based 
on a gradient method of first order. Table 1 represents the energies of the ground 
states for some light even-even nuclei. We can see that the results given by the 
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FMD model with the SLy4 parametrisation are comparable to full HF calculations, 
and that both are close to the experimental data. A similar quality of results is ob- 



E GS (MeV) 


FMD 


HF° 


Exp Y 


4 He 


-26.156 


-26.7 


-28.296 


6 Li 


-30.225 


-32.478 


-31.99 


8 Bc 


-44.996 


-45.439 


-56.5 


12 C 


-77.4 


-90.578 


-92.163 


16Q 


-127.9 


-128.485 


-127.62 


40 Ca 


-344.98 


-344.22 


-342.056 



Table 1. FMD ground state energy with the SLy4 interaction for some even-even nuclei in com- 
parison with results from the Hartree-Fock model, and experimental data. 

tained for mean square radii. A discrepancy is however observed for 12 C where FMD 
underestimates the binding energy. This can be tentatively understood observing 
that the FMD result has a pronounced a structure; similar states is experimentally 
observed as excited states. Another interesting observable is the GMR frequency, 
which is a stringent test of the compressibility properties of nuclear matter. To 
excite a Giant Resonances, at a given time an exciting operator is applied to the 
system according to : 

H = -^Nuclei + XOS(t) (7) 

where O = r 2 for the GMR, which characterizes the isotropic compression of the sys- 
tem at the initial time. In calculating the GMR for 12 C with the SLy5 parametrisa- 
tion, we found that the frequency of the monopole is given by ujfmd = 24.017MeV, 
which is comparable to the experimental value luexp — 21.9 ± 0.3MeV extracted by 
D.H. Youngblood et alP 1 . Calculations for heavier systems and different multipo- 
larities are in progress. 
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3. FMD for thermodynamics 

We have seen that the FMD model gives a good description of light systems at 
T=0MeV, ie good ground state properties and compressibility. We will propose 
now a method to extract the thermal properties at finite temperature. 

The extraction of the phase diagram at finite temperature for a finite system 
needs the calculation of thermal averages of observable < A > a , as well as of the 
temperature T. If the system is ergodic, the thermal average is given by the time 
average < A > = linit^oo \ J A(t')dt' . To know the temperature of the system, 
this latter has to be coupled to a thermometer, if we want the microcanonical 
temperature, or to a thermal bath if we want the canonical temperature. In both 
cases the coupling must be such that the total wave function is a tensor product 
between the nuclear system and the thermometer | 'J Total) = (^Nuclear) ® I^Th), 
where the total Hamiltonian is given by : 

#Total = -^Nuclear + #Th + H\ (8) 

Here, H\ is the weak coupling between the two systems where the expectation values 
should respect the conditions < ^Nuclear >^< "Hi > and < Hrh >^><Hi >. If the 
statistical ensemble is the microcanonical ensemble (thermometer case), we should 
always respect the condition < 7YNuciear >3>< Htii >, while for the canonical 
ensemble (coupling to a thermal bath) we should have < H.Th >^< ^Nuclear >• If 
the thermal eresponse of the thermometer (or the thermal bath) < H-rh >o= f{T) 
is known analytically, it can be used to extract the common temperature of the two 
systems at equilibrium from the calculated time average < Ti-Th >o- 

In practical calculations, the nucleus is trapped in a potential, and the total 
Hamiltonian for the nuclear system can be written as : 

-^Nuclear -^Kinetic + -^Interaction + -^Trap (9) 

For the trap we have chosen a r 4 potential such as -ffTrap = W4, 2i=i where W4 is 
the trap constant. The thermometer has been taken as an harmonic oscillator where 
-^Th = 2^ + ^ m th^th^ h , with thermal response T = faj/ln ( <^>°^I^ ) • 
The coupling between the two system has been chosen as a contact interaction 
H 1 = A J2t=i & f^th ~ FA . Even if the parameters W±,\ and fiuj are chosen to opti- 
mise the thermalisation properties, the equilibration time is very long, of the order 
of 50000 fm/c. For this reason we have made calculations only for a very small 
system, namely 4 He. Applying the method explained above at different total ener- 
gies, we have extracted the temperature as a function of the excitation energy per 
nucleon, the so-called caloric curve. Figures Q] represent three calculations for a 4 He 
nucleus using an hybrid ensemble intermediate between the microcanonical and the 
canonical ensemble, with < Tlrh >~< ^Nuclear >• The caloric curves obtained with 
the SLy4d and the SGII parametrisation, with and without Coulomb interaction 
are presented. It is important to remark that the excitation energy axis contains 
the constraining potential and as such it is not a direct measure of the deposited 
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energy. At low excitation energies the trend is not clear, and doubts may be shed 
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Fig. 1. Caloric curves for 4 He with W4 = O.OlMeV.fm -4 . The black, red, and blue dots represent 
respectively the calculations with the SLy4d interaction, SLy4d + Coulomb interaction, and the 
SGII interaction. 



on the ergodicity of the dynamics. At higher energy the curves show the expected 
trend for a transition from an excited fermion system to a system of non-interacting 
classical particles. This can be further confirmed looking at the energy distribution, 
which is bimodal in the transition region. Comparing the three caloric curves, we 
can see that that the choice of the nuclear interaction has small effects, and the 
same is true for the influence of Coulomb. Calculations with heavier systems are in 
progress. 
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